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SUMMARY , 

One of the early objectives of space probes sent to Mars and Venus -will be 
to determine the characteristics of the atmospheres on those planets for scien- 
tific purposes. Characteristics of interest include profiles of gas density, 
temperature, and pressure above the planet surface, and chemical composition of 
the atmosphere. Because of practical limitations, it is doubtful that the first 
probes will convey a very complete picture of the atmosphere, and a good initial 
objective will be to define the above properties broadly, although certainly the 
more that can be learned from an early probe the better. 

So far, the entry probe has been considered merely a vehicle for delivering 
instruments to planets, and the interactions of the probe with the atmosphere 
have been considered problems to be overcome if the mission is to be accomplished. 
The present paper is based on the premise that these interactions are a fertile 
source of information about the atmosphere. The development of this premise 
suggests some experiments which are simple (and therefore it is hoped reliable), 
lightweight, and productive of valuable information on the planets' atmosphereo. 

It is found that the density, pressure, and RT product profiles in the 
atmosphere can be obtained from measurements of axial acceleration and the angle 
of entry into the atmosphere. For a certain class of entry body configurations , 
measurements of the static flight stability can be used to give information on 
the mean ratio of specific heats of the atmospheric gases. Measuring the ambient 
temperature so as to obtain the gas constant R from its product with tempera- 
ture, RT, and hence determine the mean molecular weight of the atmospheric gases 
is discussed. Luminosity detectors that view the excited high-temperature gases 
in the probe shock layer for evidence of composition are discussed as well as the 
possibility of determining the kind of gases from measurements of stagnation 
temperatures . 

An important characteristic of this kind of probe experiment is that the 
definition of atmospheric properties begins as high in the atmosphere as any 
detectable effect on the probe occurs - in general, around 300,000 feet altitude. 


INTRODUCTION 


One of the early objectives of the exploration of Mars and Venus by space 
probes will almost certainly be to determine the composition, pressure, density, 
and temperature profiles of their surrounding atmospheres. Not only is this 
information of scientific interest for what it can indicate about the nature of 
the planet, its evolutionary state, and the possibilities of supporting life but 
it is also essential for the proper design of subsequent, more advanced entry 
vehicles. 


The constraints under which this exploration must be accomplished are as 
follows. The initial instrument payload weights will be small because of booster 
limitations. The communications problem is difficult because of the large dis- 
tance to earth and the limited power available on the vehicle. The reliability 
requirement is severe since the instruments and communications equipment must 
remain in operable condition from 10 to 40 weeks during interplanetary transit, 
and then function perfectly for a brief period during which the mechanical load- 
ing and heating may be severe. These constraints suggest that the best possibil- 
ity of success lies in the choice of simple instruments and passive systems. 
Elegant and complicated equipment may give promise of defining the atmospheric 
environment precisely and in great detail, but is incompatible with small size 
and weight, simple communication, and long-term reliability. 

In the design of the entry probe vehicle, numerous problems arising from the 
interaction of the vehicle with the atmosphere must be treated. These include 
the capsule flight stability, the aerodynamic heating, both radiative and convec- 
tive, and the determination of the flight trajectory. The designer, preoccupied 
with the solution of these problems and hampered by uncertainty in the nature of 
the atmosphere to be entered, is inclined to view them as obstacles to the ful- 
fillment of the mission. However, one may also view the vehicle interactions 
with the atmosphere as a primary source of data about the nature of the atmos- 
phere. As a simple example previously considered for application in the earth's 
atmosphere, the deceleration of the vehicle combined with a known drag coeffi- 
cient and speed will give the atmospheric density. There are other similar possi- 
bilities which, if compatible with the requirements of simplicity, reliability, 
and^Iow weight (and, in the author f s opinion, they tend to be), should perhaps be 
seriously contemplated for the early probe missions. 

In this leport, a number ol possible ways in which interactions of a probe 
with the atmosphere can lead to knowledge of the atmosphere are described. 

Methods are included for experimentally determining atmospheric density, pressure, 
and RT product profiles, mean ratio of specific heats, ambient temperature, 
mean ^molecular weight of the atmospheric gases, and percentage of carbon dioxide 
in mixtures of carbon dioxide and nitrogen. The emphasis is on concepts and 
principles. No attempt is made to evaluate the relative merit of the experiments 
described, although it is certain that some of the measurements described will be 
found to be simpler, or more useful, or more feasible than others. Without fur- 
ther detailed engineering and scientific study, however, it is not possible to 
say conclusively what the relative merits are. The principal purpose of the pres- 
ent paper is to put these techniques for planetary atmospheric research before 
the interested reader. 
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NOTATION 


a 


a x 


A 

B 


CD 


Cm 

Cnia 

Cin ai 

Cp 

d 

D 

f 

F 


h 


hs 


I 

m 

M 


speed of sound 

acceleration in flight direction 

reference area for gas "dynamic coefficients 

CpA po 

ballistic constant, ~ (3 sin 0 

drag force 
drag coefficient, (p/gJpj^A 

re storing moment 

static -stability coefficient, (l/2)p V^Ad 

dCm 

static -stability derivative, 

static -stability derivative at a = 0 

P " Poo 

pressure coefficient, ^ 1 / 2 ) p 

reference length for gas -dynamic coefficients 
acceleration profile function (appendix A) 
frequency of pitching oscillation 

frequency constant characteristic of probe (appendix A) 

gravitational acceleration within the atmosphere of the planet under 
exploration 

altitude above planet surface 
stagnation-point enthalpy of the flow 

moment of inertia about transverse axis through center of gravity 
vehicle mass 

V 

flight Mach number, — 


M.W. molecular weight of atmospheric gas mixture 
p static pressure 

R gas constant for atmospheric gas mixture, ^ 


R u universal gas constant, R X M.W. 

t time 

T gas temperature 

V flight velocity 

W vehicle weight on earth 

a angle of attack, angle between vehicle axis and flight direction 

3 reciprocal of atmospheric scale height, ^ 

7 ratio of specific heats of atmospheric gas mixture 

p gas density 

ct radius of gyration of probe about a transverse axis through the center of 

gravity 

6 flight -path angle measured relative to horizontal 

£ gas -dynamic damping coefficient (appendix A) 

Subscripts 

00 ambient conditions in undisturbed atmosphere 

e earth reference values 

E conditions at entry into the atmosphere 

s at the flow stagnation point 

0 at the planet surface 

u upper portion of the atmosphere 

1 lower portion of the atmosphere 
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EXPERIMENTS BASED ON ENTRY BODY DYNAMICS 


Atmospheric Density 


The density of a planet's atmosphere may, in principle, be determined for 
any instant of flight from a measurement of the deceleration. The governing 
equation is simply Newton's second law of motion with the defining equation for 
the drag coefficient representing the drag force. 


n 1 „ 2 . dV 
C D 2 P v A = " m dit 

(la) 

2dV/dt 

(lb) 

p V 2 CpA/m 


The mass and frontal area of the entry vehicle are known prior to flight, and the 
velocity must be determined. The velocity at entry is obtainable from deep -space 
trajectory data provided by tracking and machine computation. The velocity after 
entry may be determined from integration of the accelerometer output. If the 
measurement is restricted to altitudes where continuum flow occurs, and this 
is expected to be the greatest part of the atmosphere, the drag coefficient may 
be considered known. Experimental facilities on earth are capable of defining 
the drag coefficient for wide ranges of speed, Reynolds number and attitude, and 
for various gases. Presently available experimental data (ref. 1) indicate that 
the drag coefficients of several types of entry bodies are insensitive to the 
composition of the atmospheres through which they fly. Hence the accuracy of the 
determination should not be seriously limited by present uncertainty in the com- 
position of the planets' atmospheres. However, if the composition could be 
determined better from other results of a probe flight or from other experiments, 
the drag coefficient for the proper gaseous environment would then be measured in 
ground facilities . 

The effect of angle -of -attack oscillations on the measurement could be 
significant. Since the accelerometer would most conveniently measure the axial 
force rather than the drag, and since the axial force tends in general to dimin- 
ish with increasing angle of attack, a proper interpretation of the data would 
require angle -of -attack sensors in the probe. These are desired for another 
reason as well, and will be discussed in a later section. 

The principal limitation on accuracy of the measurement of density would 
probably be in determining the velocity from accelerometer data. The accelerom- 
eter should be able to sense accelerations over a tenfold range with an accuracy 
of the order of 1 percent of the reading. This accuracy is required because the 
output of the accelerometer must be integrated to obtain the velocity. If meas- 
urements are continued to the time when V/Ve = 0.1, if the error in measured 
acceleration is systematic and constant at 1 percent of the reading, and if there 
is also a 1-percent error in Ve which adds to the other error, then the error 
in density will increase as the entry progresses and will reach a maximum of 
approximately J+0 percent. The following table shows some accuracies calculated 
for these hypothetical conditions, as a function of the instantaneous flight 
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velocity ratio. 


v/v E 


Percent error in p 


1.0 

.8 

•5 

.2 

.1 


3 

4 

7 

20 

ko 


The altitudes at which these velocity ratios are obtained will be considered in a 
later section. It is evident that an accurate accelerometer is needed to main- 
tain accuracy in density throughout the entry. 

An alternative way of measuring atmospheric density at high vehicle veloci- 
ties is by means of the pitot pressure. The pitot pressure coefficient is 
defined by 

n _ Ps Poo , . 

°Ps - (l/2)p^V 2 

where p s is the pressure recorded at a stagnation-point orifice and p is the 
ambient pressure in the free stream. At high speeds, p^ may be neglected, 
having values about one -thousandth of p g . Hence the pitot -pressure equation may 
be simplified and solved for the density to obtain 

2p s 

p = Cp s v 2 (3) 


Thus the ambient density is defined by a pressure measurement, given the velocity 
and the coefficient Cp g . (it will be noted; however, that the accelerometer is 
still required to provide the velocity history. ) 

The pitot -pressure coefficient, like the drag coefficients of blunt bodies, 
is insensitive to gas composition at high speeds. The limiting value given by 
Newtonian theory, 2, represents the condition of total streamwise momentum loss 
at the bow wave and is a good approximation at high speeds for both ideal and 
real gases. In the following table, theoretical values of Cp are given for 
various gases, characterized by different values of the ratio of specific heats, 
and for various speeds. 
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c Ps 

i. 4 o 

5 

1.8l 

i. 4 o 

00 

I.85 

1.67 

00 

I.76 

1.00 

CO 

2.00 

Real air 

> 10 

1-9 
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It is apparent that a value of 1-90, selected without knowledge of the gas com- 
position, would be correct within 5 percent; and this would be subject to itera- 
tive improvement from straightforward real -gas computations, given a composition 
of the atmosphere. 

It may be remarked that a measurement of the pitot pressure is essentially a 
measurement of the local drag in the neighborhood of the stagnation point; and 
that it does not differ in principle from the measurement of drag as a way of 
determining ambient density. Pressure measurements at other locations of the 
orifice - not at a stagnation point - can also be used to obtain free-stream den- 
sity, provided the local pressure coefficient is known. For example, pressure 
orifices could be used on cones with attached bow waves. The orifice location 
should in any case be chosen so that the measured pressures will be large com- 
pared to free-stream pressure, to minimize the error introduced by taking the 
difference, p - n , in the numerator of equation (2). 

^ j. - 00 


Atmospheric Pressure and RT Product 


From a measured history of density as a function of time, the local static 
pressure may he deduced with the aid of the barometric equation 


1 dT 1 dp ^ 

T dh p dh “ RT V ' 

where g/RT = (3, the reciprocal of the scale height. If the atmosphere Is nearly 
isothermal, so that 


dT « d£ 

T p 

(percent change in temperature with altitude is small compared to percent change, 
in density), equation (4) reduces to 


dp gp 

dh RT 


which may be put in the form 


RT 


gpV sin 0 
dp/dt 


(5) 

(6) 


It follows that the local static pressure in the atmosphere may be obtained from 


P 


pRT 


gp 2 V sin 9 
dp/dt ” 


(7) 


Measured quantities are the density, rate of increase of density, instantaneous 
flight velocity, and flight -path angle. 
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Another, somewhat more direct, way to obtain the static pressure from the 
measurements of density is to integrate the basic barometric equation 

dp = -gpdh ( 8 ) 

to obtain 

nh 

P = / (-gp)cLh (9) 

'Ooo 

Since 

dh = -(V sin 0 )dt, 
nt 

P = / (gpV sin 0 )dt ( 10 ) 

^ o 

This amounts to summing up the weight of gas above the given altitude. The dif- 
ferential form, equation (7), works on a local basis, first obtaining RT from 
local measurements of p and dp/dt. The integral form, equation (lO), has one 
apparent advantage - it does not require the assumption of an isothermal atmos- 
phere. However, this assumption is normally not a serious one, so that both 
methods may be of use in analyzing a set of probe data. 

From an accuracy analysis, based on the assumed performance of the acceler- 
ometer cited earlier, it is concluded that pressure can be obtained more accu- 
rately than the density. This follows since in equation (lO) pV is known more 
accurately than p (see eq. (lb)). Also, the accuracies obtainable by the two 
evaluation techniques (eqs. (7) and (lO)) are not the same. The following table 
gives some calculated values. 


Percent error in p 


v/v E 

Equation (7) 

Equation (lO) 

1.0 

2.1 

2.1 

.8 

1.4 

2.5 

.5 

1.5 

4 

.2 

19 

10 


r+3^9 

20 

.1 

1 -55 


Equation (7) is the more accurate for early stages of the entry, but becomes very 
bad at the low velocities. The favorable error levels for equation (7) at higher 
velocities were calculated with the assumption that there were no random noise - 
type errors in the acceleration measurement - only the systematic errors 
described earlier. Noise would, of course, degrade the definition of the deriva- 
tive dp/dt and would tend to favor the use of equation (lO). 


8 



The RT product defined hy equation (6) is also basically useful informa- 
tion. The determination of either R or T by any technique then permits the 
other quantity to be known also. For examp le, if T is measured, R may be 
obtained, and, in turn, the mean molecular weight 

M.W. = TT (U) 


of the atmospheric gases may be calculated. However , its accuracy is limited by 
the accuracy of R, which is approximately the same as the accuracy of p 
obtained by the integration method. A discussion of possible measurement tech- 
niques for obtaining atmospheric temperature will be given later. 


In the above analysis no error in determining pressure was attributed to 
definition of sin 0 . In both equations (7) and (10), however, the flight -path 
angle 9 is needed. It is also required to convert p(t) and p(t) to altitude 
■profiles p(h) and p(h), since time after entry is not a fundamentally important 
variable and what we really wish to know is the dependence of p, p, etc., on 
altitude. The relation between time and altitude is as follows: 


h 


f t impact 


( V sin 0 )dt 


,J t 


Altitude may be calculated from the flight -path angle 


( 12 ) 


9 (t ) and the velocity V(t). 


One possible way of satisfying the need to know sin 9 in equations (7 )j ^ 
(10), and (12) is to guide the entry probe into a nearly vertical entry, 6 ~ 90 , 
so that sin 9 ~ 1. The actual entry angle may then depart as far as ±lb from 
vertical and still introduce only 5 -percent error if sin 9 is assumed equal to 
1.0. This error might be considered acceptable in equations (7)> (10), and. ( 12 ), 
in which case, a guidance system capable of bringing the probe into a vertical 
entry within ±l8° would be required. 


A more general and perhaps more satisfactory solution to this problem would 
utilize the simplest possible angle sensing instrumentation in the probe. One 
system of this kind relies on the assumption that, regardless of entry angle, 0 
approaches 90° near the end of flight. This tends to be true whenever the probe 
slows to a very low speed before impact, because the gravitational acceleration ^ 
then becomes the primary cause of motion. If it can be assumed that ©final 90 
and that the probe remains gas -dynamically stable until impact, then a simple 
gyroscope in the probe to indicate angle change from entry to impact would give 
the entry angle. The gyroscope would be alined with the body axis just. prior to 
entry, and then allowed to remain fixed in inertial space. The angle history of 
the probe axis relative to the gyro axis might then be like. that shown in 
sketch (a). This system would define 0-g and also the attitude angle and pitch- 
ing frequency of the vehicle. Further study would be required to determine 
whether such a system would adequately fill the need for flight -path -angle data. 


Another way to determine 0, that in some respects is very attractive, is^to 
carry on board the entry vehicle a radio altimeter. This direct measurement ox 




Some estimates, which have come to the writer's attention, however, indicate that 
this instrumentation is excessively heavy and bulky. It is beyond the scope of 
■the present note to judge the merit of its inclusion. 

V in ^o PC l nt in thiS connection should be made. It was noted that entries 
near 0-90 minimize the accuracy with which 0 must be determined. The oppo- 
si e is also true flat entries near the overshoot boundary require the most 
precise measurement of 0, and also can lead to "skipping" trajectories in which 
/dt is first negative, then positive, and then negative again. The determina- 
tion of atmospheric pressure by equation (10) would be considerably complicated 
y this type of trajectory. It is felt that steep entries free of skipping are 
preferable for early flights, for reasons of simplicity and accuracy. 


Studies of Gas Composition From Static 
Flight Stability Measurements 


The gas -dynamic moments responsible for static flight stability, being the 
residual moments from a balance of positive and negative moments applied over the 
vehicle, tend to indicate with high sensitivity small changes in the distribution 
of pressures. Thus the stability may be affected to a considerably greater 
extent than the drag by changes in gas composition. One type of configuration 
for which this is true is shown in figure 1. The data, taken from reference 1 
show that the stability derivative increases progressively as the ratio of spe- 
cific heats increases. Note that the flight stability of this body in helium and 
argon correlates along a single curve, in spite of the large differences in molec- 
ular weight and gas constant. This indicates that the ratio of specific heats is 
the gas property primarily responsible for the changes in stability. 


From a measurement of the Mach number and the stability derivative of this 
configuration (or one having similar behavior), it should be possible to estimate 
the mean ratio of specific heats in the planetary atmosphere. Since the Mach 
number l'^ )0 (t) is up to this point not determined, it is necessary to discuss in 
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some detail how this analysis would he made. It must he done as a simultaneous 
solution using the experimental curves of figure 1 and a second known relation 
between Mach number and 7. Consider that a value of Cj^ has been obtained. 
This measurement has a locus on figure 1 which is a level line, since Mach number 
is unknown. Along the level line, a definite variation of 7 with M holds, 
which is just a cross plot of the data in figure 1, as shown in sketch (b) by the 
line labeled = constant. A second 

relation between M and 7 is provided by 
the definition of M, 

M = — == ( 14 ) 

V7 RT 

Here, RT and V are known by the proce- 
dures of the previous section. Hence, we 
may write, for any time of the flight 

7M 2 = known constant 

which is also plotted schematically in 
sketch (h). The intersection of the two 
curves is the solution point and gives 
the desired values of 7 and M at time 
t. The analysis is repeated for other 
times. It would be expected, of course, 
that for an atmosphere of fixed composition, the same value of 7 would be 

obtained at all times . 

The stability derivative is obtained from measurements of the frequency of 
the attitude oscillation provided by the attitude -indicating gyroscope discussed 
earlier or its equivalent. The defining equation is 



2rtt 


/( “Cm a )(l/ 2 )pV 2 Ad 


( 15 ) 


which, when solved for 
tion, becomes 


Cma as a function of acceleration in the flight direc- 


Cma = ^ 



"D 


ff 

ax 


(16) 


Thus, the stability derivative is given by known constants, the drag coefficient, 
the pitching frequency, and the axial acceleration. 

Although the curves in figure 1 have, for simplicity, been labeled with the 
value of 7 , ratio of specific heats in the free stream, it is well known that 
7 varies with position in high-temperature flow fields as a result of dissocia- 
tion, vibration, and other chemical and thermodynamic processes. Hence, it is 
probably more precise to consider that every gas mixture has, at a given ambient 
pressure, a characteristic curve in the coordinates of figure 1 . Thus, all gases 
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of ?oo = 1,11 would - not produce the same stability curve , because of differences 
in dissociation temperature, products of dissociation, etc* This may be regarded 
as an advantage in that it provides, in principle, a method for distinguishing 
between gases of the same For this reason, obtaining the trace of stability 

as a i unction of Mach number over the widest possible range of Mach numbers 
should be valuable. 

It is worth noting that attached flow over blunt-nosed slender bodies is 
reasonably well understood theoretically (refs. 2 and 3). The effects of gas com- 
position on stability calculated with this theory roughly agree with measurements, 
although the stability derivative cannot be obtained theoretically with the pre- 
cision necessary for these experiments at the present time. The theoretical 
understanding would, however, be of value. For example, it might be necessary to 
consider possible effects of nonequilibrium chemistry on the stability for gases, 
such as carbon dioxide, which when dissociated are very slow to recombine. 

Other problems would require experimental support, notably cases involving 
separated flow. Also, experiments should be relied on to supply the dependence 
of stability on Mach number for comparison with the flight data. Additional 
ground -facility experiments could define these curves better for the gases which 
have been considered and for other gas mixtures. 

It seems probable that other configurations will be found that are sensitive 
in static stability to gas composition. It would be particularly desirable to 
find a body which has this property but is not subject to flow separation effects, 
such as can occur on the flare of the present configuration. It may be that 
bodies for which the afterbody pressures are dominated by the blast -wave produc- 
ing influence of a blunt nose will, in general, show some degree of j 
sens itivity. 


OTHER MEASUREMENTS 
Atmospheric Temperature 


Aside from its basic value as a property of the planet's atmosphere, the 
atmospheric temperature has the potentiality of leading through equation (ll) to 
the determination of the mean molecular weight of the atmospheric gases as noted 
earlier. Unfortunately, the free -stream static ten^perature of the atmosphere has 
essentially no gas -dynamic effect other than to determine the free-stream Mach 
number. Thus, Mach number and velocity measurements used in conjunction with 
equation (l4) would merely supply a redundant value of RT and would not assist 
in defining R or T separately. 

Near the end of the flight, when the speed becomes subsonic, the gas 
temperatures are not much affected by fluid-dynamic effects and the atmosphere 
temperature may be measured directly to permit the determination of R. (This 
would require precautions to avoid measurement errors due to the presence of the 
hot heat shield.) At high speeds, however, this is not possible because of com- 
pression and frictional heating of the gas around the probe. It may be suffi- 
cient to defer atmospheric temperature measurements to the low-speed period and 
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to determine the mean molecular weight only at low altitudes. It would then be 
reasonable to infer temperature at higher altitudes by assuming constant R, 
given RT(h). Unfortunately, however, the accuracy in RT measurements deterio- 
rates at low speeds, so that such a procedure might result in a serious loss of 
accuracy in determining the molecular weight. 


Two techniques have occurred to the author for determining atmospheric 
temperature in the free stream at the higher speeds and altitudes. While they do 
not appear extremely attractive, the end result of determining the mean molecular 
weight accurately and as a function of altitude may be judged to be of such great 
value that these or equivalent techniques should be considered. The first is to 
conduct on-board experiments with a gas sample as schematically indicated in 
sketch (c). The gas sample is drawn through a cooler into a test chamber where 
(a) its temperature is measured, and (b) 

its sound speed is determined from the Microphone Temperature sensor 

frequency of reverberation of a single 
pulse introduced by the microphone, which 
serves both to excite and to measure the 
vibrations. The gas constant and molecu- 
lar weight are then computed from 

R = a 2 /7T = Ru/M.W. (17) 


V 




Test chamber 




-Cooler 


t 


Sample 

Sketch (c) 


inlet 


A serious difficulty with this approach is 
the possibility of contamination of the 
gas sample by heat -shield ablation prod- 
ucts, which, it is believed, would be dif- 
ficult to avoid. Another source of error 
is the possibility of dissociated atmos- 
pheric gases being frozen in the dissoci- 
ated condition while being cooled. A 
further problem is the small time avail- 
able for readying the sample, particularly 
if the measurement is to be made at sev- 
eral altitudes, since the entire entry 
period is less than half a minute (see 
Trajectory Considerations). These consid- 
erations make this technique appear some- 
thing less than promising. 

A second method is to release a low- 
density subprobe designed to decelerate 
to low speed at a high altitude, to meas- 
ure the static temperature directly and 
telemeter it to the entry probe for recording and later relay. For example, a 
subprobe weighing 3 pounds with a frontal area of 1 square foot and a drag coef- 
ficient of 1.5 (m/CpA = l/l 6 slug/ft 2 ) would decelerate from an entry velocity. of 
26,000 ft/sec to I 3 OO ft/sec at an altitude of 220,000 feet in the Jet Propulsion 
Laboratory Mars atmosphere "E. " It would then take several minutes to drift down 
to the planet. It might have a configuration like that shown in sketch (d). A 



Sketch ( d) 
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telemetry range of pO miles would be required from the subprobe to the probe for 
steep entry angles , and much of the transmission would occur after the main probe 
landed. Measurement error due to the elevated temperature of the subprobe heat 
shield could be avoided by discarding the heat shield after the subprobe slowed 
down. This approach, while somewhat complicated, perhaps offers a chance for 
determining atmospheric temperature and molecular weight profiles at high 
altitude . 


Shock-Layer Spectroscopy 


Because of its high speed, the probe entering the atmosphere generates a 
shock layer of high -temperature gas that is luminous at speeds above about 
-^-3^000 ft/sec* The radiation is normally considered as a source of heat input to 
the capsule, but it can also provide data on the atmosphere composition. The 
shock compression ahead of a blunt nose performs essentially the same function as 
the spark, arc, or flame used in spectroscopy - it excites the constituent gases 
and causes them to emit characteristic radiation, which, under proper conditions, 
can be interpreted to indicate which gases are present. 

The shock layer would be contaminated with gases evolving from the heat- 
shield ablation. The ablation vapor radiation has a characteristic spectrum, 
which depends primarily on the heat-shield materials but may depend somewhat on 
the atmospheric composition. It has been determined from tests at the Ames 
Research Center, at speeds between 20,000 and 30,000 ft/sec, that the predominant 
ablation product radiation for many materials is in the infrared region of the 
spectrum, reflecting their lower temperature, while that from the atmospheric 
gases peaks in the ultraviolet. Hence, it appears that contaminant radiations 
may be avoided by judicious selection of materials and of wavelengths for the 
detectors . 

The effect of varying the percent of C0 2 in nitrogen and carbon dioxide 
mixtures of about the proportions believed to occur on Mars and Venus is illus- 
trated in figure 2. These data are from reference 1, and were obtained in the 
Ames pilot hypersonic free -flight facility at a mean speed of 21,500 ft/sec and 
an ambient pressure of 0.02 atmosphere. Although data were obtained over the 
spectral range from 2,000 to 10,000 A , the points shown in this figure are from 
a single ultraviolet detector centered at 3,660 A with a band width to half sen- 
sitivity of 380 A. A very sensitive response to small amounts of carbon dioxide 
is indicated. Sensitivity to mixture composition is also shown by comparison of 
these data with data for air, given in this figure just to the left ofs£he origin. 
The radiation detected here is believed to originate from cyanogen and from ion- 
ized nitrogen molecules. 

The best wavelengths for detectors to perform spectroscopic analysis of a 
planet’s atmosphere would certainly depend on the atmosphere composition assumed. 
For N 2 and C0 2 mixtures , it would appear that one detector should respond to the 
ultraviolet band as shown in figure 2; it should perhaps be made more selective 
to separate the cyanogen and the nitrogen molecule radiations. Minor constitu- 
ents such as water vapor, oxygen, or argon could also be investigated by shock- 
layer spectroscopy. It would appear desirable also to employ one broad-band 
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detector to measure total luminosity for comparison with expected total luminos- 
ity. This would serve as an over -all check on the correctness of the composi- 
tion, since the presence of gases other than those presently postulated might be 
missed in any narrow -band detector, but would be expected to affect the level ^of 
the over -all radiation. The output of the broad -band detector would necessarily 
be corrected for radiation from ablation products on the basis of laboratory 
tests such as those referred to above. 

The results obtained from this kind of instrumentation would be strengthened 
by the fact that they produce a line trace on a map of luminous intensity as a 
function of speed and density, rather than an isolated point of data. The trace 
would be compared with a reference map obtained from laboratory tests in atmos- 
pheres of assumed composition at the flight velocities and densities. Agreement 
with the reference map at all times of flight would be rather convincing evidence 
of the correctness of the assumed composition. Furthermore, the necessity for 
agreement at all points of the map should assist in ruling out incorrect 
interpretations of the data. 


Stagnation Temperature Effects 


The stagnation temperature of a gas mixture is a sensitive function of the 
gas composition at a given ambient pressure and total enthalpy, or equivalently, 
at a given flight velocity and altitude. Temperature is shown in figure 3 as a 
function of enthalpy for pure carbon dioxide and pure nitrogen for pressures of 
0.01 atmosphere and 1 atmosphere. The data for this figure were taken from ref- 
erences 4 and 5* The corresponding curves of stagnation temperature as a func- 
tion of flight velocity are shown in figure 4. The large variations in kinetic 
temperature with gas composition are noteworthy. At velocities between 10,000 
and 20,000 ft/sec, stagnation temperature differences as large as 4500 R occur 
between these two pure gases. 

Unfortunately, this gas temperature is difficult to measure. Direct meas- 
urement with thermocouples must be ruled out for Up, and for C0p in the early 
parts of the flight, because the temperatures are high enough to vaporize the 
sensor. The radiation measurements described above constitute a temperature 
measurement, in a sense, but the temperature becomes known only when the composi- 
tion is established. What is desired here is an independent, more direct 
determination. 

Ordinarily, a convective heat -transfer measurement is a way of indirectly 
obtaining the driving temperature potential. However, in reference 6, it is con- 
tended that heat transfer in dissociating gas systems is driven by the total 
enthalpy and not by the total temperature. The writer would argue that this can- 
not be true in the case of elevated wall temperatures , because in the limit where 
wall temperature equals stagnation temperature, no heat transfer due to either 
conduction or diffusion would be expected. Hence, the wall temperature for zero 
heat transfer is a strong function of the gas stagnation temperature rather than 
of stagnation enthalpy. For wall temperatures which are high but still below 
material temperature limits, important differences may be found in heat transfer 
for a given total enthalpy but for various total temperatures - that is, for 
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various gas compositions. This problem requires analysis by methods similar to 
those of reference 6, and hence requires extensive numerical computations on elec 
tronic computers. This was considered beyond the scope of the present paper, but 
would certainly be of interest as a possible method of determining gas composi- 
tion from convective heat transfer measured as a function of flight velocity. 


Trajectory Considerations 


For measurements of the kind considered above, the characteristics of entry 
trajectories are of interest and importance since they indicate the altitude 
range over which measurements will be made, the time period within which all data 
will be collected, the frequency and convergence of pitching oscillations, the 
accelerations to which the vehicle and the instruments will be subjected, etc. 

The author undertook to calculate some of these details of the trajectories using 
the method of Allen and Eggers, reference 7* Although this method is more approx 
imate than some others that are available, it shows directly the effects of the 
principal variables and is quick and easy to apply. It is also acceptably accu- 
rate for steep-angled trajectories, which are not close to an overshoot boundary, 
that is, trajectories which are nearly straight-line paths. This is shown by fig 
ure 5 for some nearly vertical entries into Mars and Venus where the approximate 
calculations are compared with numerical integrations of the complete motion equa 
tions. Certainly the lack of knowledge of the atmospheres greatly overshadows 
any disagreement shown in figure 5. 

The atmospheres used in these estimates were some of those selected by the 
Jet Propulsion Laboratory, NASA, for planetary studies, based on information 
given in references 8, 9> a nd 10- They are shown on logarithmic coordinates in 
figure 6 with the straight-line fairings used in the present calculations to 
approximate them. 

Equations used to obtain the following results are given in the appendix. 
They were either taken from reference 7 or can be easily derived from those given 
in the reference. 

Velocity -altitude histories are shown in figure 7 for both Mars and Venus. 
Values of the ballistic coefficient m/Cj^A, ranging from l/l6 to 6 slugs/ft 2 , 
entry angles of 30 and 90 , and various atmospheric structures are included on 
this figure. For Mars, the following conclusion emerges (contingent on the cor- 
rectness of these atmospheres): Unless measurements are made during the probe 

entry, a major part of the sensible atmosphere of the planet will not be investi- 
gated. In the worst case, the probe penetrates to an altitude of 26,000 feet 
before it has slowed to 1000 ft/sec, so that low-speed measurements, made from a 
parachute, for example, would be confined to the near vicinity of the surface. 

From the gas -dynamic standpoint, as can be seen from this figure, the atmosphere 
begins at about 300,000 feet. It will certainly be of interest and importance to 
explore this upper region. An alternative to using the vehicle entry for probe 
purposes which still permits relatively high altitude measurements is to make the 
vehicle very light so that it decelerates to low speed high in the atmosphere. 

The curve for m/Cj)A of l/b indicates this method of approach. 
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For the Venus atmosphere, the limitation described above is not so severe. 
With reasonable vehicle parameters, the atmosphere could be •explored by low-speed 
(parachute) techniques below an altitude of 170,000 feet. The region from there 
to 300,000 feet would again, however, be lost to this approach. 

The time period of entry, as noted in the appendix, is directly proportional 
to the scale height of the atmosphere, l/p, and inversely proportional to 
Vg sin 0. The entry times are shown in very general form in figure 8 as a func- 
tion of the scale height for entry angles of 3O 0 and 90°. This plot is made for 
Ve of 26,000 ft/sec, but can easily be modified for any other entry speed. The 
times shown are between speeds of 0.95 Ve and 0.05 Ve* Extension of the velocity 
interval on either end greatly increases the time interval, but this is somewhat 
misleading since the principal deceleration occurs in the periods shown. For 
vertical entry at scale heights as small as 20,000 feet, the entry times are very 
short - in the neighborhood of 10 seconds . Such short times pose a requirement 
for the instrumentation systems that they be fast -acting. Of course, most of the 
instruments considered above can be made to have submillisecond response times, 
so perhaps the main requirement imposed by the short entry time is that all the 
data be tape recorded for later telemetry to earth at a more deliberate speed. 

The tape recording is also required for another reason - the well-known communi- 
cations blackout phenomenon due to flow ionization which can be expected to occur 
at high speeds . 1 


The short entry times are naturally accompanied by high peak decelerations. 
These are shown for Mars entry speed in the lower part of figure 8. These values 
indicate that the instruments must be of good strength -to -weight ratio to remain 
undamaged during entry. 

Calculations of the frequency and damping indicate frequencies from 1 to 
10 cycles per second, and highly convergent oscillations (see equations in appen- 
dix). Bodies considered here were 2 feet in diameter. The frequency profile 
function, which shows the relative variation of pitching frequency with velocity 
during entry and which is universal within the limitations of this analysis, is 
shown in figure 9* and typical oscillatory amplitude variations are shown in fig- 
ure 10. These calculations indicate that entries with small pitching rate at 
initial angles as great as A5^ should pose no particular problems for measurement 
of stability. 


CONCLUDING REMARKS 


It appears that considerable information about the atmospheres of the 
planets could be obtained from measurements of entry probe interactions with the 
atmospheres. This information can be divided into two classes: one pertaining 

1 The high atmospheric temperatures - up to 7^0° K - presently believed to 
occur low in the atmosphere of Venus suggest that the tape-recorded data should 
be transmitted as rapidly as possible to the ,T bus n vehicle from which the probe 
was launched and stored there for transmission to earth. The author has been 
informed that the designs of electronic systems employing transistors to operate 
at elevated temperatures as high as 7^0° K is difficult if not impossible. 
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to the pressure, density , and RT product profiles in the gas above the planet 
surface, and the other pertaining to composition. The latter could not be 
defined completely or exactly, it is felt, particularly if the atmosphere is a 
complex mixture of gases similar to that on the earth. However, the principal 
constituents might be identified, the presence of some minor constituents might 
be investigated, and the mean molecular weight determined, which would certainly 
be a large gain over present knowledge. 


It is evident that much work would be required in a number of areas to 
develop a probe of this type. Ground research would be required to aid in select- 
ing important design variaDles. Instrument developments would be required, and 
additional error analyses would have to be made. This paper is far from complete 
in respect to these matters. It is hoped, however, that it will stimulate inter- 
est in the use of gas -dynamic effects for investigating the atmospheres of the 
planets. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Jan. 23 , 1963 
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APPENDIX A 


EQUATIONS USED FOR ESTIMATING ENTRY TRAJECTORIES 


A basic equation of reference 7 which relates velocity to altitude (density) 


is 


where 


V 2 
% = e 


B - 

P 


(Al) 


B = 


c dA Pc 


m p sin 0 


E 


and 


Equation (Al) in the form 


p -ph 
p = — = e 

Po 


2 , V E 

p = b Zn T 


is used to eliminate density in favor of velocity in the frequency equation, etc. 

For the lower segment of an atmosphere represented by two intersecting 
exponential density variations 


p = e 


-Puh 


h > hi , p Q = p Q 


‘U 


-p z h 

p = e 0 


h < hi , .p Q = p 


P = b7 ln ¥ + p i h < hl 


(A2) 


where Vi and p are conditions at hi. 


The time period of entry is 


t = 


1 r dh 

sin Q E J V 


which, by use of the preceding relations, can be put in the form 
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t(v) 


i r VE/ v d(v R /v) 

3V E sin e E Jv E /Vi ^E/V) 


(A3) 


Equation (A3) gives the time elapsed in passing from velocity V± to velocity 
Note that it is independent of m/CpA and p Q . With the integral evaluated for 
Vi/\ E - 0-95 and V/V E = 0.05, equation (A 3 ) is plotted in figure 8 . 


V. 


The acceleration is given by 


C D (l/2)pV 2 A 


which can be reduced to the form 

a x = PV E 2 sin 0 E 



again independent of m/CpA 


and p o . 



The function 


(£)*■* 


V 


(A4) 


which describes the acceleration profile during entry, is plotted in figure 11 . 
At any altitude , the instantaneous pitching frequency is 

£ _ _^I_ / ( ) ( 1 / 2 ) pv* Ad 

2ir V I 


which, in terms of the above variables , can be written 


f 


= FV E 


/Dp sin 0 E 
v d 


where D d(V/V e ) is the acceleration profile function, d is the reference 
length for the pitching moment^ and 


p = r 7i f 

2n: v Cp \aj 

contains the aerodynamic and inertial characteristics of the model. The pitching 

frequency profile is just the square root of the acceleration profile and has been 
shown in figure 9 . " ' 
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For estimating the pitching amplitude variations, the pitching motion 
equation from reference 11 is employed 


a 

a E 


„ , J 

A e\ t 

VdN ( Cma\ in(V E /V) 1 

W 0 

L\oy \ c D y pa sin e E _ 


1 / 2 


where J Q is the zero-order Bessel function and a E is the. angle of attack at 
entry. For present estimates , the aerodynamic damping function 

i = C D - C]^ + (Ci^ + Cn^Kd/a ) 2 


was taken ^ ^1 to zero. This function is strongly dependent on configuration, 
taking on both positive and negative values, but does not greatly influence the 
angle history until low velocities are reached, near the end of flight (see 
refs. 11 and 12). 


After the second positive peak of the Bessel function, which occurs at 
J (x) = 0.300, x = 7-0, V/V E ~ 0«99> the following approximation for the amplitude 

of pitching holds 


a 

a E 


0.300 



t/ 2C D, 


ln(V E /V 2 ) 

Zn(V E /V) 


1/4 


v 2 being the exact value of velocity at the second positive peak of the Bessel 
function. 
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30 x 10 



Light emitted near 3700 A wavelength per cubic centimeter of gas in the shock layer as 
function of amount of C0 2 in N 2 -C0 2 mixtures. 









Figure 4. - Stagnation temperature as a function of flight velocity for pure 

and pure C0 2 . 










= 0.00238 siugs/ ft 3 



Figure 6.- Model atmospheres used for trajectory calculations. 
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Figure 7 • " Altitude -velocity relations 
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Figure 7 • - Concluded 






















